#brokenStick function
#ARGS:
#	n - number of broken stick values
#RETURN:
#	broken stick distribution vector
brokenStick <- function(n){
	valVec <- 1/c(1:n)
	resultList <- lapply(c(1:n),function(x){
				sum(valVec[x:n])/n
			})
	return(unlist(resultList))
	
} 